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SYMBOLS 


[A] stiffness matrices for discretized finite-element region 

[a],[b] upper triangular matrices 

[A q ],[B q ] element stiffness and mass matrix, respectively, for element q 

a width of rectangular element 




acoustic pressure function for impedance tube 

respective value of pressure eigenfunction and mean-flow profile at node Z 
one-dimensional shape functions 

height and width, respectively, of impedance tube 

free-space wave number 

dimensional axial propagation constant 

dimensionless axial propagation constant 

node number 

mean-flow profile 

centerline Mach number 

cross-mode order 


n integer 1 or 2 denoting particle velocity or particle displacement, 

respectively 

NY, NZ total number of elements in y- and z-direction, respectively 

N two-dimensional shape function 

Z 

SPL sound pressure level 

P acoustic pressure 


iii 



t 


time 


x,y,z Cartesian coordinates 

x ,y ,z Cartesian coordinate values at node I 

acoust ^ c admittances of wall lining 
acoustic admittance 

dimensionless local coordinate system for element 
eigenvalue for cross-mode order, m 
global vector of unknowns 

vector of unknown nodal values for element q 
angular frequency 
Subscripts: 

I,J,m integers 

Superscript: 

tr vector transpose 

Mathematical notation: 

[ ] matrix 

{ } vector 

det | | determinate of matrix 


n,5 

X 

m 

(♦i 

i+ q } 

a) 
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INTRODUCTION 


The reduction of high intensity noise radiating from jet engine aircraft has 
been the subject of continuous investigation for the past two decades. At the pres- 
ent time, mathematical models for predicting the noise radiated from these engines 
are well developed (ref. 1). However, these prediction models cannot be used with 
confidence to predict noise levels radiating from a jet engine in a realistic flow 
environment. Before this confidence can be achieved, it is necessary to obtain accu- 
rate determinations of the physical parameters, which are input to the prediction 
models. The two critical input parameters are the source input and the admittance of 
the sound-absorbing material. This paper addresses the determination of the latter 
of these two important input parameters. 

Experiments with sound-absorbing materials in the presence of flow indicate 
significant differences in acoustic performance as the mean-flow Mach number is var- 
ied (ref. 2). This necessitates that sound-absorbing properties of acoustic materi- 
als be evaluated experimentally in flow. To account for the effects of the mean 
flow, a test sample of the acoustic-lining material is normally installed in a 
grazing-flow impedance tube with a controlled aeroacoustic environment. In such a 
grazing-flow impedance tube, usually rectangular in cross section, the acoustic mate- 
rial is aligned so that the sound and mean-flow graze over its surface in a manner 
simulating an engine environment. In principle, the acoustic admittance can be 
determined by direct measurements of normal particle velocity and pressure since 
acoustic admittance is the ratio of normal particle velocity to pressure at the wall 
surface. However, measurements of normal velocity in the presence of grazing flow 
are not reliable. Instead, an indirect measurement approach (ref. 3) has been used 
which makes use of changes in the wave field caused by the introduction of the test 
sample in the hard-walled impedance tube. These measurements are then used to calcu- 
late the acoustic admittance. 

Analytical methods for directly determining acoustic admittances from pressure 
measurements taken in a grazing-flow impedance tube have progressed slowly. 

Armstrong, Beckemeyer, and Olsen (ref. 2) developed a method that is restricted to 
thin boundary-layer flows with shear in one cross-sectional direction of the tube. 
This method used an asymptotic expansion to derive a boundary condition that is 
applicable at the boundaries of the uniform region in which an exact solution is 
possible. Mungur and Gladwell (ref. 3) developed a method that is applicable to more 
general boundary layers, including cases where the boundary layer spans the entire 
tube. This method is based upon a Runge-Kutta integration across the boundary layer 
and is restricted to grazing flows with shear in a single cross-sectional direction. 
Watson (ref. 4) developed a finite-element method for calculating acoustic admit- 
tances from acoustic pressure data in a grazing-flow impedance tube with one- 
dimensional shear . 

The methods presented in references 2, 3, and 4 are restricted to infinitely 
long test samples with mean shear in a single cross-sectional direction of an imped- 
ance tube. However, grazing-flow duct facilities commonly employ ducts with a rela- 
tively small rectangular cross section in which the grazing flow possesses shear in 
both cross-sectional directions of the impedance tube, such as the flow impedance 
test laboratory at the Langley Research Center. Thus, the present effort was moti- 



vated by the need to account for the more realistic flow environment with shear in 
both cross-sectional directions of laboratory flow-impedance tubes. 

This paper describes a method for determining the admittance of an acoustic 
material in a grazing-flow impedance tube in which the grazing flow is assumed to 
possess shear in both cross-sectional directions. The analysis is restricted to 
sound fields which are dominated by a single propagating mode. The method is devel- 
oped explicitly for rectangular ducts, although the approach is applicable to other 
duct geometries as well. Experimental input data needed include an estimate of the 
axial propagation constant, the two-dimensional mean-flow profile, and acoustic pres- 
sures parallel to the treated wall of the impedance tube. From these experimentally 
determined characteristics of the aeroacoustic field, the admittance can be calcu- 
lated. In this paper, the unknown admittance value is obtained by solving an eigen- 
value problem. This eigenvalue problem results from the application of the finite- 
element method to the partial differential equation and boundary conditions governing 
the acoustic field. 

The paper is divided into four sections. The first section presents the gov- 
erning partial differential equation and boundary conditions. In the second section, 
a finite-element scheme is applied to obtain the solution to this problem. The 
finite-element method leads to an eigenvalue problem which requires a special scheme 
to obtain its solution. This eigenvalue solution scheme is developed in the third 
section of this paper. Finally, the paper closes with a discussion of some 
applications. 


MATHEMATICAL FORMULATION 

Figure 1 depicts the cross-sectional geometry of the grazing-flow impedance-tube 
test section and its coordinate system. All four walls of the tube are lined with 
acoustic treatment and are assumed infinite in extent. The mathematical development 
assumes that the mean flow in the test section is axial and fully developed, so that 
flow gradients exist only in the y- and z-direction. 

The solution of the wave equation for harmonic time dependence in the presence 
of shear can be written in the form: 


-i (K x x-ojt) 

P(x,y,z,t) = F(y,z)e 


( 1 ) 


where is the axial propagation constant. The following elliptic partial differ- 

ential equation for the pressure function F(y,z) is obtained (ref. 5): 


E(y,z) = (l - MK x )[F yy + F zz ] + 2K x [M y F y + M Z F Z ] 
+ K 2 [(l - MK X ) 3 - K^(l - MK X )]F = 0 


( 2 ) 
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where 


“x 


K 


In equation (2), the function M = M(y,z) is the mean flow profile, K x is the 
dimensionless axial propagation constant in the impedance tube, and K is the free- 
space wave number (w/c). 

The physical boundary conditions associated with equation (2) require either 
continuity of particle displacement or continuity of normal particle velocity along 
the four walls of the impedance tube (ref. 1 ). These boundary conditions can be 
expressed in the form: 


F 2 (y,0) = -iK&Jl - M(y,0)K x ] n F(y,0) (3) 
F z (y,L) = iK$ 2 [ 1 - M(y,L)K x ] n F(y,L) (4) 
F y (0,z) - -iK$ 3 [l - M(0,z)K x ] n F(0,z ) (5) 
F y (H,z) = iK 6 4 [l - M(H,z)K x ] n F(H,z) ( 6 ) 


where 3-], & 2 ’ ^ 3 ' an< ^ 3 4 are t ^ ie specific acoustic admittances of the acousti- 

cally lined walls. (See fig. 1.) Continuity of acoustic particle displacement at 
the lined walls are obtained by setting the integer n to 2 , whereas continuity of 
acoustic particle velocity is obtained by setting n to 1 . It is instructive to 
note that realistic mean-flow profiles satisfy the condition of no slip at the bound- 
aries rendering a zero mean flow there. Thus, for realistic mean flows, both forms 
of the boundary condition coalesce to the same expression at each wall. However, 
since some unrealistic mean velocity profiles are used in this paper in order to 
check with former works, n will be carried throughout the analysis as a parameter. 

Equations (2) through ( 6 ) constitute a boundary value problem for the acoustic 
pressure function F(y,z). The axial propagation constant k , free-space wave num- 
ber K, mean-flow profile M(y,z), boundary condition order n, and wall-lining 
admittances 6 ^ , 3 2 / an< ^ 83 are assumed specified and are treated as known. The 

homogeneity of equations ( 2 ) through ( 6 ) then allows the determination of the func- 
tion F(y,z) and unknown acoustic admittance 3 . 

4 

Upon obtaining the value of the necessary parameters, the boundary value problem 
posed by equations (2) through ( 6 ) can be solved to obtain the function F(y,z) and 
unknown admittance 3 4 > However, equation (2) is a partial differential equation 
with variable coefficients and its solution can be put in terms of known functions 
only for some special cases of the two-dimensional sheared flow, M(y,z). However, 
these cases cannot generally be achieved in the laboratory and are not useful for 
general application. Thus, the admittance and pressure eigenfunction F(y,z) must 
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be determined numerically for mean flows of practical interest. A numerical method 
for finding these solutions is described in the following section. 


NUMERICAL SOLUTION TECHNIQUE 

Two principal decisions in the adoption of a numerical technique were the use of 
a Galerkin finite-element method and use of linear shape functions to approximate the 
pressure function and mean-flow profile. The discretization mesh is shown in fig- 
ure 2. The cross-sectional area of the impedance tube has been divided into NY 
and NZ elements in the y- and z -direction, respectively. A typical element (I,J) 
has length (z I+1 - z I ) and (y^^ - Yj) in the z- and y-direction, respectively. In 
order to generalize and also simplify the formulation, it is convenient to define a 
local coordinate system for each element. This local coordinate system, shown in 
figure 3, also facilitates the integration which is required to obtain the element 
equation. 

In order to develop the element equations for a typical element (I,J), the pres- 
sure function F(n,5) and mean-flow profile M(n*S) are expressed in a form which 
ensures continuity of these functions both within the element and also along lines 
which are common between any two of the elements. In a typical element, these func- 
tions are approximated by linear interpolation functions of the form 


F(n,5) = N £ <n,S) F £ 
M(n,5) = N £ (n,£) m £ 


U * i to 4) 


a = i to 4) 


(7) 

( 8 ) 


in which repeated subscripts are to be summed over, F £ and M £ are the values of 
the functions F(ti,£) and M(n,£) at node l of the element, and 


N 1 = f 1 (n,b) f 1 (5, a) (9) 

N 2 = f^rwb) f 2 U,a) (10) 

N 3 = f 2 (n ' b) f 2 (£ ' a) (11) 

N 4 - f 2 (n,b) ( 12 ) 


f 1 (n,b) = i - £- 


4 


(13) 



f 2 (n,b) 


n. 

b 


(14) 


In applying the Galerkin version of the method of weighted residuals, the resid- 
ual error E(n,5) in each element is multiplied by the basis function N £ (n,£) and 
integrated over the element to obtain 



E(n,5) N A (n,5) dn d£ 


This integral expresses the desired averaging of the residual error but it does not 
admit the influence of the wall boundary conditions. Integration by parts is now 
incorporated to reduce the order of the second derivative terms in this integration 
and to incorporate the boundary conditions: 



-(/.‘ft 


E(n,S) N £ (n,?) dn d? 


9F Hk + IE _ K 
3*0 3n 3£ 35 x \3ri Urj 


+ 3? H 


a - m 


3M 3F x , \ , 3M 3 F M 


+ K 2 ^ - K^)mFN £ + K 2 (k£ - 1 ) FN £ 

- 3K 2 K^MFN £ + K 2 K^M 3 FN £ J dr) dgj 

+ J j[l - K x M(b,?)] N £ (b,z) 

- [l - K X M(0,5)] |f(0,5) N £ (0,5)j d£ 
+ J" j^1 - K x M(n,a)] -|jp(n,a) N £ (n,a) 


- [i - K x M(n,0)] ||<n,0) N £ (n,0)j 


dn 


(15) 
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Boundary conditions are now incorporated at the element level by replacing the normal 
derivative of the function F along the outer boundaries of the discretized region 
by its value from the boundary conditions. For example, T^*(b,5) along the boundary 
y * H is replaced by iK3 4 [l - M(b,£)K x ] n F(b,£). This procedure is standard prac- 
tice in finite-element methodology since it avoids the difficult task of having to 
choose basis functions which satisfy the boundary conditions. Substituting equa- 
tions (7) through (14) into equation (15) gives 


j J" E(n,C) N^(n,5) <35 <3n = [A q ]{<|> q } + 6 4 [ B q ] {<|> q } 06) 

in which [A q ] and [B q ] are the stiffness and mass matrix for the element, g^ is 
the unknown admittance value, and {<(> q } is a vector containing the unknown values of 
the function F(t}, 5) at the four nodes of the element. 

Assembly of a global matrix equation representing the discretized region 
(fig. 2) from a set of finite -element matrix equations is a basic procedure in the 
finite-element technique. Appropriate shifting of rows and columns are all that is 
required to add the local element matrix directly into the global matrix (ref. 6). 

The result is a set of linear matrix equations which is expressible in the form: 


[A] {<|)}= e 4 !B]{*} (17) 

The global matrices [A] and [B] are both square matrices whose order is 

(NY + 1 ) (NZ + 1) and {$} is a vector consisting of the values of F(r),£) at the 

various nodes of the system. 

Equation (17) has the form of a generalized eigenvalue problem. Here, the 
unknown admittance 3^ is the eigenvalue and is determined by the determinant 
condition 

det I [A] - 3 [B] I « 0 (18) 

4 

Equation (18) is the basic equation which must be solved to obtain the unknown 

admittance 3 „ 

4. 

Some further comments concerning equation (18) should be made in order to avoid 
some confusion. First, there are many values of 3 4 which satisfy equation (18) and 
this appears to be at odds with the physical problem for which only a single value of 
the admittance 3 4 exist. Secondly, equation (18) must be solved with an eigenvalue 
extraction algorithm on a digital computer. In this paper, the LZ algorithm (ref. 7) 
is employed on a Control Data 6600 computer system to exact the solution to equa- 
tion (18). The order of the coefficient matrices [A] and [B] was restricted to 
100 when the algorithm was employed on this machine. This restriction is due to the 
core limitation of the central memory of the computer and is equivalent to a maximum 
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of nine elements in both the y- and z-directions of the impedance tube. Since 10 
elements per wavelength are needed to resolve the acoustic wave, this would restrict 
the analysis to only plane waves, a situation which rarely exists physically. 

In the following section, a determinate search routine is developed which takes 
advantage of the character of the coefficient matrices in equation (18) and increases 
the maximum number of points in the y- and z-direction of the impedance tube by more 
than an order of magnitude. It is also shown how to extract the correct value of the 
unknown admittance B ^ from the many values which are obtained from the solution to 
equation (18). 


DETERMINATE SEARCH TECHNIQUE 

Generalized eigenvalue problems for large systems such as equation (17) are 
solved by numerical schemes that are either direct or iterative. With the improve- 
ment in digital computer hardware, direct methods have proved to be more versatile 
and reliable and are' used in this paper. In the direct approach, both [A] and 
[B] are transformed to upper triangular form and the solution to the transformed 
system is obtained directly. Fortunately, in the finite-element application, the 
equations are amendable to direct solution techniques that take advantage of the 
special character of [A] and [B] . The general form for the coefficient matrices 
[A] and [B] are as follows: 


a 1,1 

a 1 , 2 



a 2,1 

a 2, 2 

a 2 , 3 



iuJ 

tsj 

a 3 , 3 

a 3 , 4 


[A] = 


(19) 


a NY+1 ,NY a NY+1,NY+1 


0 0 


0 0 0 


[B] - 


0 0 


0 


( 20 ) 


0 


b NY+1 , NY+1 
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While the unknown vector {<)>} is partitioned in the form (eq. (17) ): 


{*} 


tr 


tv V V •**' *NY' ^NY+1 1 


(21) 


In equations (19) through (21) the submatrices a^ j and ^ N y+1 NY+1 are s< 3 uare 
matrices whose orders are (NZ + 1), and are* column vectors of the same order. 

Further, the submatrices a^j and b NY+1 NY+1 are tr ^- d ^ a 9 onal * ’Thus, [A] is a 
block tridiagonal matrix and [B] has a single nonzero tridiagonal block in its last 
row and column. 

Because of the character of the coefficient matrices [A] and [B] , equa- 
tion (17) is fairly easy to solve. There are three steps that are executed 
recursively: 

1 . Use fully pivoted row operations to reduce [A] to block upper triangular 
form, executing the same row operations on the matrix [B] . The resulting general- 
ized eigenvalue problem will be of the form: 


[a] {<#>} = e 4 [i]{«j>} 


( 22 ) 


a 1,1 a 1 , 2 

a 2, 2 a 2, 3 



a 3, 3 a 3,4 


a NY,NY 3 NY, NY+1 

a NY+1 ,NY+1 


(23) 
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[i] , 


0 

0 0 
0 


NY, NY+1 


J NY+1 ,NY+1 


(24) 


2. The eigenvalue is now extracted from the last equation of this reduced 

system (eq. (22)): 


[ a NY+1 , NY+1 -I ^NY+1 ^ ^4 t b NY+1 ,NY+1 ^ ^NY+1 ^ 

Equation (25) is solved by using the LZ algorithm (ref. 7) to obtain and 

{ ^ NY+ 1 1 * Tbe rema ^ n; *- n 9 unknowns, {<t> NY K ^NY -1 ^ ' •••• {‘I’.,} are obtained by back 


(25) 


[ S NY , NY -I ^ ^NY ^ ^4 b NY,NY+1 S NY,NY+1 ^ ^NY+1 ^ 


(26) 


11 - sv ' OT - v ' •••' 


2) 


(27) 


3. Matrices [A] and [a] are written to mass storage in rectangular form and_ 
the nonzero blocks in [B] and [b] are stored into the zero blocks of [A] and [a] 


[A] = 


0 

a 1,1 

a 1 , 2 

a 2, 1 

a 2, 2 

a 2, 3 

a 3 , 2 

• 

a 3 , 3 

• 

a 3 , 4 

• 

• 

a NY,NY-1 

« 

a NY,NY 

• 

a NY,NY+1 

a NY+1 , NY 

a NY+1 , NY+1 

b NY+1 , NY+1 


(28) 
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b NY , NY+ 1 

a 1,1 

a 1 , 2 

a 2, 1 

a 2, 2 

a 2, 3 

a 3, 1 

a 3, 3 

a 3,4 

• 

. 

. 

9 

. 

. 

a NY,NY-1 

a NY,NY 

a NY,NY+1 

a NY,NY 

a NY+1 , NY+1 

u 

NY+1 , NY+1 


( 29 ) 


The advantage of the technique employed herein is 

1 • The [B] matrix need not be stored* 

2* Only enough space for three columns of the [A] matrix is needed in mass 
storage* 

3. The number of points in the z-direction of the impedance tube which the 

scheme can now handle is approximately 100 instead of 10, and the number 
of points in the y-direction of the tube is unlimited* 

It should be noted that the solution to equation (25) gives the admittance g^ and 
unknown vector However, the vector { 1 ^ contains values of the 

function F(y,z) at nodes along the upper wall of the discretized region (fig. 2)* 

To obtain values of the function F(y,z) at points away from the upper wall, equa- 
tions (26) and (27) must be solved* The thrust of this paper is to determine the 
unknown admittance g^; therefore, only equation (25) need be solved. 

Some further discussion on the solution to equation (25) is helpful* As men- 
tioned earlier, there are many values of g^ which satisfy this equation. In fact, 
there are (NZ + 1) values of g^ which satisfy this equation, and there will be in 
general a different eigenfunction {^ny+I^ for each g^* This result is not 
surprising since equation (2) is elliptic in character and the boundary condition 
along the upper wall, y = H, was not specified (i.e., the boundary condition was not 
specified along a closed contour). In short, the discrete set of wall admittances 
obtained from the solution to equation (25) simply reflect the fact that a unique 
solution to equations (2) through (6) cannot be obtained unless additional boundary 
information is given along the upper wall* The solution to equation (25) gives all 
possible combinations of the upper wall eigenfunctions {^y+l} anc * wa ^ admittance 
values which can satisfy the boundary condition in equation (6). Suppose a measure- 
ment of the upper wall eigenfunction {^^y+il made the impedance tube. Then 
this measured upper wall eigenfunction will correspond to one of the solutions for 
{^NY+1 1 obtained from solving equation (25). The admittance g^ corresponding to 
this eigenfunction is the correct one ror the problem. This approach will be used in 
this paper. 
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RESULTS AND DISCUSSION 


In order to verify the numerical technique, results from the method are compared 
both with exact solutions that can be obtained for a uniform mean-flow profile and 
with results of reference 4 for cases involving shear in the y-direction only. Next, 
results are presented for some cases in which the flow possesses shear in both the y- 
and z-direction of the tube. Results are restricted to a square duct for which 
H « L = 1 foot, K = 1 (radians per foot), with 50 elements in the y- and 
z-directions of the tube (NY = NZ = 50). Further, continuity of particle velocity is 
employed as the boundary condition and the bottom and side walls of the tube are con- 
sidered rigid (B-j = $2 = ^3 = * 0,0i). The choice of this particular combination 

of parameters was felt to be sufficient to obtain confidence in the numerical 
technique • 


Uniform Flow 

The fundamental equations governing the acoustic pressure eigenfunction 
(eqs. (2) through (6)) can be solved analytically under the assumption of uniform 
mean flow. This expression for the eigenfunction is expressed in the following form: 


— / » __ , \ mirz 

F(y,z) = F (y) cos — — 
m L 


(m — 0, 1 , 2, ) 


(30) 


cos X y 
„ , N nr 

F (v) = 

m * cos X H 
m 


(31) 


where the eigenvalue X is related to the axial wave number by the equation: 

m 


X 


2 

m 



i)k| - 2M^ + 


2 2 
_ m tt 

K 2 L 2 , 


(32) 


Equations (30) through (32) are easily derived by finding the general solution to 
equations (2) through (5) for a constant M by separation of variables. Note that 
equation (30) satisfies equations (3) through (5) for rigid boundaries (i.e., 

3<| = $2 = ^3 ~ + A generalization of this solution to include nonrigid 

boundaries is also possible (ref. 1). Equation (30) is incomplete by virtue of the 
fact that it does not satisfy the upper wall boundary condition. (See eq. (6).) The 
substitution of equation (30) into equation (6) gives the following equation: 


*4 = 


-X w tan (X H) 
m m 


iK( 1 - MK ) 
x 


n 


1 1 


(33) 



Generally, H, L, K, M, m, n, and 3^ are specified. The objective is then to 
find the discrete set of axial propagation constants which satisfy equation (33). A 
single infinity of these axial propagation constants can be found. In this paper, we 
want to consider the inverse problem. Thus, we wish to specify H, L, K, M, n, 
and k^; then it is a simple matter to calculate the corresponding values of 3^ 
from equation (33). There will be a single infinity of these values, one for each 
value of m, so that 


B 4 - (5 


0 ' *1 


5 2' *•*} 


(34) 




m 


-X m tanU H) 
m m 

iK( 1 - MK x ) n 


(35) 


Note that for the problem considered here, a unique value of 3 4 is not defined 
until a prescription is given for singling out a single value for the integer m. 

The prescription employed is to specify the functional form of the eigenfunction 
F(y,z) along the upper wall. Evaluating equation (30) along the upper walls gives 


F ( H, z ) = cos (m = 0, 1, 2, ..., °°) 

L 


(36) 


Thus, a specification of F(H,z) is equivalent to specifying a single value of m. 

The major thrust of the numerical method employed here is embodied in the solu- 
tion to equation (25). The solution to equation (25) gives the discrete set of 
admittance values specified in equation (34). Further, {^Y+1^ so ^ ut ^ on f° r 

F(H,z). If the numerical method has been applied properly, each of the upper wall 
eigenfunctions Hny+i} computed from equation (25) will be of the form cos irnrz/L, 
just as in equation (36). In the rest of this section, these ideas will be confirmed 
with numerical results. 

Consider first the solution to equation (35) for K x = 2.97819i and zero flow 

(M = 0) . In this case, considerable insight into the integrity of the numerical 

scheme can be obtained since the solution to equation (35) gives 3 = 0 + Oi for 

m < 2, and 3 purely imaginary for m > 2. Although this particular case has no 
m 

relevance to physical reality, it does test the ability of the numerical method to 
extract both zero and multiple roots simultaneously. Results are presented in 
table I, starting at m = 0 and terminating at m = 5. Excellent comparison between 
the exact admittance values and those extracted by the finite-element method are 
obtained. The wall pressure eigenfunction F(H,z) is plotted against the exact 
value for m = 0, m = 1, and m = 2 in figures 4(a), (b), and (c), respectively. 

The wall pressure eigenfunctions computed from this analysis are undistinguishable 
from the exact values. 


1 2 



Table II presents results for K x = 0,5 + 0.5i and zero flow. In this case, 
the admittance 8 posesses both real and imaginary parts. As seen in the table, 
agreement between m the exact admittance value and that computed from the finite- 
element method are good. The wall pressure eigenfunctions for m * 0, m = 1, and 
m = 2 are plotted in figures 5(a), (b), and (c), respectively. Excellent comparison 
is again obtained. 

A sample calculation for uniform flow, in which M = 0.2 and K x = 0.5 + 0.5i 
is given in table III. Numerical values of the wall admittance are observed to be in 
good agreement with the exact values given in the table. The wall pressure eigen- 
function, F(H,z), for m = 0, m * 1, and m = 2 were also compared with their 
exact value. Although the plots are not shown for the sake of brevity, comparisons 
were consistent with those in figure 5. 

The final case involving uniform flow is for a Mach 0.5 flow (M = 0.5), with 
K x = 0.5+ 0.5i. Results for the predicted admittance are tabulated in table IV, 
and the accuracy is good when compared with the exact values. Wall pressure eigen- 
functions for m - 0, in = 1, and m = 2 for this 0.5 Mach number case were plotted 
and even though results are not shown, trends were identical to those of figure 4. 


One-Dimensional Shear 

8 m 3M 

When there is one-dimensional shear so that ^ = 0 and * 0, an exact 

expression for the admittance is apparently not possible. However, because there is 
no variation of the mean-flow profile in the z-direction of the tube, numerical re- 
sults can be compared with those of a one-dimensional analysis presented by Watson 
(ref. 4). The results which are presented from reference 4 use 1000 points in the 
y-direction. Further, sample calculations are restricted to a constant gradient 
sheared flow for which 


M (y) 



(o < r < in 
(!«*<■)] 


(37) 


where M Q is the centerline Mach number in the tube. 

Admittance values for a centerline Mach number of 0.2 (M q = 0.2) and 
K x = 0.5 + 0.5i are tabulated in table V. Quantitatively, the eigenvalues computed 
from the two-dimensional analysis presented in this paper compare well with those 
obtained by the one-dimensional analysis presented in reference 4. The wall eigen- 
functions F(H,z) corresponding to the admittances tabulated in table V have also 
been computed. Because the mean-flow profile is independent of the coordinate z, ' 
the eigenfunction F(y,z) is still defined by equation (30). The wall pressure 
eigenfunction along the upper wall will be of the form, cos imtz/L, just as in the 
uniform flow case. However, the transverse eigenfunction F m (y) and wa ll admittance 
values 84 are no longer given by the analytical expressions defined in equa- 
tions (31) and 33). In fact, F (y) must be determined numerically. The wall pres- 
sure eigenfunction F(H,z) computed from this analysis are compared with the exact 
values for m = 0 , m = 1 , and m = 2 in figures 6 (a), (b), and (c), respectively. 
Values computed from this analysis are indistinguishable from the exact values. Such 
good comparison between the modes gives further credence to the analysis employed in 
this paper. 
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The final sample calculation is for M = 0*5 and K x = 0.5 + 0.5i. Admittance 
calculations are given in table VI. Note that for this reasonably high centerline 
Mach number, the admittance values computed from this analysis compare well with 
those computed from reference 4. The wall pressure eigenfunction, F(H,z), for 
m = 0, m = 1, and m = 2 were consistent with those in figure 6. 


Two-Dimensional Shear 

Results are now presented for a flow containing gradients in both cross- 
sectional directions of an impedance tube. In this final example, a demonstration of 
how data taken from a two-dimensional sheared flow experiment are used in conjunction 
with the analysis presented here to determine the admittance of a test specimen is 
given. The experimental data were obtained in the flow impedance test laboratory at 
the Langley Research Center. 

A sketch of the grazing flow impedance tube is shown in figure 7. The tube is 
2 inches square, so that H = L = 2 inches, and the data were obtained at 2500 hertz. 
Note that the origin of the coordinate system is chosen in the upper left-hand corner 
of the tube, so that the test specimen is located on the upper wall when referenced 
to this origin. Data obtained from the axially transversing microphone establish the 
attenuation and phase rate characteristics along the length of the test specimen. 
Attenuation and phase rate data are translated directly into an axial propagation 
constant ( ref • 2 ) . 

A plot of sound pressure level (SPL) in decibels as a function of axial distance 
at 2500 Hz and M = 0 is shown in figure 8. The leading and trailing edge of the 
test specimen is denoted by an asterisk on the X-axis. It is of interest to note 
that the slope of the curve in figure 8 at axial position x gives the attenuation 
rate at the position. Further, the curve is linear in regions which are dominated by 
a single propagating mode. The standing wave patterns in the vicinity of the leading 
and trailing edge of the specimen indicate the presence of reflections and higher 
order mode contamination. The analysis presented in this paper is not applicable in 
these regions because of the presence of multiple modes in the acoustic field. How- 
ever, the curve is clearly linear for values of x between 27 and 35 centimeters. 

It is this linear portion of the curve that was used to obtain the attenuation rate. 
Although the curve is not shown, phase rate data show trends similar to that in fig- 
ure 8 with the linear portion of the curve covering a much wider range of values 
of x. 


The mean-flow profile in the tube was obtained by aerodynamic measurements, at 
the 350 points in the cross section. The node points for the finite elements were 
chosen at points in the domain where mean flow measurements were taken (i.e., NZ = 9 
and NY = 34). Figure 9 gives a graphical representation of the flow profile mea- 
sured in the tube as a function of the coordinate y, at z = 1, 1/2, and 1/4 inch. 
Note that the functional form of the mean-flow profile changes at each of the three 
z stations so that this profile is not independent of the coordinate z as in the 
previous section. Further, the centerline Mach number of this profile was determined 
to be approximately 0.3 (M q = 0.3). 

A sample calculation is shown in table VII for no flow, uniform flow, a flow 
with one-dimensional shear, and the measured flow. The admittance given in the table 
corresponds to the admittance of the mode which is closest to a plane wave. The flow 
with one-dimensional shear was constructed by using the profile at z »■ L/2 in fig- 


14 



ure 9 (i.e., M(y,z) = G,(y)), whereas the uniform flow value was chosen as the 

•H \ 


average value of G^(y) 


(i.e., H-*j[ 


G^(y) dy j. Results show that uniform flow 


calculations agree better with results for one-dimensional shear than the no-flow 
calculation. Further, the uniform flow profile with the particle displacement 
continuity (n « 2) as the boundary condition, agrees better with the one-dimensional 
shear results than the particle velocity boundary condition (n « 1 ) . Note that the 
uniform flow and one -dimensional sheared flow admittance value differ only moderately 
from the zero flow value. In contrast the imaginary part of the admittance calcu- 
lated by using the measured profile, differs significantly from that calculated from 
the profile with only one-dimensional shear. Such results suggest that for flows 
over acoustically lined panels, the flow has a significant effect on the acoustic 
properties of the liner. The reader is reminded however, that this result would be 
invalidated if the measured value of the wall pressure eigenfunction F(H,z) did not 
coincide with that calculated from the two-dimensional sheared flow analysis. 
Unfortunately, measurements of the upper wall eigenfunction were not available and 
could not be compared with values computed from this analysis. 


CONCLUDING REMARKS 

A method has been developed for calculating the acoustic admittance of a test 
specimen located within the wall of a grazing flow impedance tube with flow gradients 
in both cross-sectional directions. The method, developed and tested for rectangular 
tubes, may be extended to other geometries as well. In this approach, the unknown 
admittance value is assumed constant and is obtained by solving an eigenvalue prob- 
lem. This eigenvalue problem results from the application of the finite-element 
method to the partial differential equation and boundary conditions governing the 
acoustic field. A major effort has been devoted toward the solution to this eigen- 
value problem, which requires a special scheme to obtain the eigenvalues. 

Admittance values determined from the method were compared with exact solutions 
obtained for a constant mean-flow profile and with results of NASA Technical Paper 
2310 for cases involving shear in a single cross-sectional direction. Excellent 
comparisons were obtained giving credibility to the scheme. Results have been given 
for the first time in which the flow is represented realistically with gradients in 
both cross-sectional directions. Limited results obtained for a test specimen 
installed in the wall of the flow impedance test laboratory at the Langley Research 
Center suggest that grazing flows can have a significant impact on the acoustic 
admittance. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
August 23, 1985 
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TABLE I.- ADMITTANCE FOR M = 0 
AND K x = 2.97819i 



Exact admittance. 

Admittance from 

m 

Pm 

finite-element method 

0 

0.00 + O.OOi 

0.00 + O.OOi 

1 

0.00 + O.OOi 

0.00 + O.OOi 

2 

0.00 - 5.44i 

0.00 - 5.45i 

3 

0.00 - 8. 89i 

0.00 - 8.98i 

4 

0.00 - 12.171 

0.00 - 12.23i 

5 

0.00 - 1 5. 39i 

0.00 - 1 5. 52i 


TABLE II.- ADMITTANCE FOR M = 0 
AND K x = 0.5 + 0.5i 


m 

Exact admittance, 
B 

m 

Admittance from 
finite-element method 

0 

1.13 + 1.211 

1.13 + 1.211 

1 

0.09 - 2.97i 

0.09 - 2.97i 

2 

0.04 - 6. 20i 

0.04 - 6. 21 i 

3 

0.03 - 9. 37i 

0.03 - 9.40i 

4 

0.02 - 1 2 . 53i 

0.02 - 12.59i 

5 

0.02 - 1 5.68i 

0.02 - 1 5.81 i 


TABLE III.- ADMITTANCE FOR M = 0.2 
AND K x = 0.5 + 0.5i 


m 

Exact admittance, 

5 

m 

Admittance from 
finite-element method 

0 

1.20 + 0. 90i 

1.20 + 0.90i 

1 

0.49 - 3. 28i 

0.49 - 3. 28i 

2 

0.82 - 6.82i 

0.82 - 6.83i 

3 

1.18 - 10. 29i 

1.19 - 1 0. 32i 

4 

1.56 - 13.951 

1.57 - 13.83i 

5 

1.94 - 17.211 

1.95 - 1 7. 35i 






TABLE IV.- ADMITTANCE FOR M * 0.5 
AND K x - 0.5 + 0.51 


m 

Exact admittance, 

Admittance from 
finite-element method 

0 

1.25 + 0.66i 

1.25 + 0.66i 

1 

1.40 - 3.60i 

1.40 - 3.60i 

2 

2.58 - 7.46i 

2.58 - 7.47i 

3 

3.82 - 11 .26i 

3.83 - 11 .29i 

4 

5.06 - 15.041 

5.09 - 15.1 2i 

5 

6.31 - 1 8. 82i 

6.36 - 1 8. 97i 


TABLE V.- ADMITTANCE FOR M Q = 0.2, K x = 0.5 + 0.5i, 
AND CONSTANT GRADIENT ONE-DIMENSIONAL SHEARED FLOW 


m 

Admittance from 
one-dimensional 
analysis (ref. 5) 

Admittance from 
present analysis 

0 

1.19 + 

1.171 

1.12 + 

1 .20i 

1 

1 

o 

CO 

• 

o 

3.2i 

0.31 - 

3.1 7i 

2 

0 

. 

to 

1 

6.41 i 

0.27 - 

6. 43i 

3 

0.24 - 

9.59i 

0.26 - 

9.63i 

4 

0.23 - 

1 2.76i 

0.26 - 

1 2.82i 


0.22 - 

1 5.94i 

0.25 - 

16.001 


TABLE VI.- ADMITTANCE FOR M Q = 0.5, K x = 0.5 + 0.5i, 
AND CONSTANT GRADIENT ONE-DIMENSIONAL SHEARED FLOW 


m 

Admittance from 
one-dimensional 
analysis (ref. 5) 

Admittance from 
present analysis 

0 

1.27 + 

1 .22i 

1.29 

+ 

1 • 31 i 

1 

0.71 - 

3.42i 

0.74 

- 

3.44i 

2 

0.63 - 

6.71i 

0.68 

- 

6.76i 

3 

0.58 - 

9.88i 

0.64 

- 

9.96i 

4 

0.56 - 1 3. 06i 

0.63 

- 

13.1 7i 

5 

0.55 - 1 6.24i 

0.62 

- 

1 6.39i 






TABLE VII.- ADMITTANCE CALCULATED FOR 
TEST SAMPLE AT MACH NUMBER OF 0.3 


Flow profile 

Admittance 

No_f low 

K x = 1.12 - 0.25i 

0.54 + 0.61 Oi 

Uniform flow 

K = 0.89 - 0.14i 
X n = 1 
n = 2 

0.49 + 0.550i 
0.67 + 0.690i 

One -dimensional flow 
K x = 0.89 - 0.1 4i 

0.61 + 0.583i 

Measured 

K x « 0.89 - 0.1 4i 

0.60 + 0.1 50i 





(NY.1) (NY, 2) (NY,3) 
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discretization. 
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Figure 3.- Global and local coordinate system for a typical finite element* 




Figure 4.- Wall pressure eigenfunctions for 




Figure 4.- Continued. 















(a) m = 0. 

Figure 6.- Wall pressure eigenfunctions for one-dimensional sheared flow 

K = 0.5 + 0.5i and Mq = 0.2. 
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Figure 7.- Schematic of grazing flow impedance tube facility. 
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M(y,L/2)=G 3 (y) 
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